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Abstract 

We present the complete differential decay rates for the process B° — > 
J/ipK + K~ including S-wave and P-wave angular momentum states for the 
K + K~ meson pair. We examine the effect of an S-wave component on the 
determination of the CP violating phase 2/3 s . Data from the B-factories indi- 
cate that an S-wave component of about 10% may be expected in the 0(1020) 
resonance region. We find that if this contribution is ignored in the analysis 
it could cause a bias in the measured value of 2/3 s towards zero of the order of 
10%. When including the K + K~ S-wave component we observe an increase in 
the statistical error on 2j3 s by less than 15%. We also point out the possibility 
of measuring the sign of cos 2f3 s by using the interference between the K + K~ 
S-wave and P-wave amplitudes to resolve the strong phase ambiguity. We 
conclude that the S-wave component can be properly taken into account in 
the analysis. 
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1 Introduction 

The decay B° — > J/00 is a golden channel for the measurement of the B° mixing phase 
— 2f3 s which is a very sensitive probe of new physics. It has been extensively studied [TJ 
§@[t|§§[7)[8). In the decay J/00, followed by a two-body decay 0(1020) -> K+K", 
the K + K~ meson pair is in an orbital P-wave amplitude. However, in the vicinity of the 
0(1020) mass, the K + K~ system can have contributions from other partial waves. The 
same comment holds for the K + K~ system in the decay channels B° — > K + K~Kg and 
D° — > K + K~7t°. The BaBar experiment showed that in these decays the S-wave and P- 
wave contributions dominate in the mass range above threshold up to 1.1 GeV/c 2 (9 10 



In both cases there is a dominant resonant 0(1020) contribution. In addition an S-wave 
fo(980) and a non-resonant contribution are found to be necessary to describe the data. 
These results motivated us to investigate the effects of a possible S-wave contribution to 
B° — > J/0K + K~ in the 0(1020) mass region. 

In the decay B° — > J/0K + K~ the K + K~ system can only arise from a ss quark pair 
while in B° — > K + K~Kg and D° — > K + K~7r° it can have contributions from both ss and 
dd. This makes it difficult to give a quantitative estimate for the S-wave component. In 



reference 11 the S-wave K + K~ contribution under the 0(1020) peak is estimated to be 
5 — 10% for decay modes in which the K + K _ arises from an ss quark pair. In this study 
we consider an S-wave of similar magnitude and assess its impact on the determination 
of the weak mixing phase — 2(3 S . 



2 Time-dependent angular distributions in the decay 
J/0K + K~ including S-wave contributions 

We consider P- and S-wave amplitudes in the decay B° — > J/0>K + K~ where the invariant 
mass of the K + K~ meson pair is in the 0(1020) mass region and the J/0 meson decays into a 
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pair. The S-wave contribution can be non-resonant or due to the fo(980) resonance^ 
We denote decay amplitudes for the B° — ► J/0K + K~ by A = (A Q , A\\, A±, As). Here A , 
An and A± are the three P-wave amplitudes consistent with the K + K system decaying 
via the 0(1020) resonance. As is the amplitude for a possible S-wave contribution in the 
K + K~ system. The amplitudes for the conjugate decay B° — > J/0K + K~ are denoted by 
A = (Aq, A\\, A±, As), which, in the absence of direct CP violation, are related to A by 
A = A , A\\ = A\\, A± = —Aj_ and As = —As- Note that A and A\\ are CP-even 
whereas A± and A s are CP-odd. The amplitudes (Aq,A\\,A±) and the amplitude A s 
may have different dependences on the mass wk+k- of the K + K~ system. However, in 
sufficiently small bins of m K + K -, such as the narrow mass region around the 0(1020) 
resonance, the dependences of the amplitudes on mx+K- can be neglected. 

We define the total P-wave strength, A 2 P = \A \ 2 + \A\\\ 2 + \A±\ 2 , the longitudinal and 
perpendicular polarisation fractions relative to the P-wave strength Rn = \A\\\ 2 /A 2 P , and 
R± = |^4_i_| 2 A4p) an d the S-wave fraction, R s = \Ag\ 2 /(A P + |^4^ | 2 ) . The phases of these 
decay amplitudes are defined by Aj = \Aj\e tSj , where j = 0, 1 1, _L, S. As only the relative 
strong phase differences can be measured we adopt the convention 5 = 0. 

An angular analysis is required to disentangle the different CP eigenstates on a sta- 
tistical basis. The angular observables are denoted as the helicity angles Q = (9i,9k,p)- 
Here 9i is the angle between the fi + momentum and the direction opposite to the B° 
momentum in the J/ip rest frame; 9k is the angle between the K + momentum and the 
direction opposite to the B° momentum in the rest frame of the K + K~ system; (p is the 
angle between the decay planes of the J/0 — > fi + n~ and the K + K~ pair, when going from 
the positive kaon to the positive lepton with a rotation around the opposite direction of 
the B° momentum in the J/0 rest frame. 

The differential decay rate for a B° meson produced at time t — decaying as B° — > 
J/0K + K~ at proper time t is given by 

d 4 r(B s ^ j/^k+k-) ™ 

T, , n J I A OcJ2hk(t)fk(ty, (2.1) 

at a cos 9 a cos ip cup rp. 
whereas the differential decay rate for an initial B° meson is given by 

d 4 r(B s ° -> MK+K-) ™ , N „ ^ 
at a cos 9 a cos yj dip 

Each of the hk(t), hk{t) and /fe(fi) for k — 1 — 10 are defined in Table [lj In total there are 
four amplitude-squared terms for the three polarisations of the P-waves and the S-wave 
component plus six interference terms. 

The time-dependence of the ten functions hk(t) for an initial B° meson state can be 
written as: 

\A (t)\ 2 = |^ | 2 e" rst cosh ^^^^j -cos $sinh ^^^^j + sin $sin(Am s t) , (2.3) 
|A||(t)| 2 = |A||| 2 e" rs * cosh (^7^ -cos $sinh (J^^J + sin *sin(Am s t) , (2.4) 



1 The mass dependence of the fo(980) is distorted as the central value of the resonance is below 
threshold. 
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Table 1: Definition of the functions hk(i), hk(i) and fk{^ii^Ki ( f) °f Eq. 2.1 and 2.2 
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where $ = — 2/3 s , Am s , Ar s and T s denote the weak mixing phase, mass difference, decay 
width difference and average decay width of the B°-Bg system, respectively. Here we have 
assumed that each of the decay amplitudes in A is dominated by a single weak phase, 
therefore a common effective 2/3 s can be used for all CP eigenstates. The time evolution 
functions hk(t) for an initial B° meson can be obtained by reversing the sign of each term 
proportional to sin(Am s t) or cos(Am s t) in hk(t). 
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3 Measuring 2f3 s in the presence of a K + K S-wave 

In this section we investigate how the measurement of 2/3 s is affected by the presence of 
a possible K + K~ S-wave contribution. We use Monte Carlo simulated toy data based 
on the differential decay rate expressions of Section [2j We generate signal decays only 
and ignore backgrounds underneath the B|? mass peak as well as all detector effects. The 
inclusion of these effects does not alter the qualitative results of this study. 



We assume a tagging efficiency e tag = 56% and a wrong tag probability u t , 



<><:) 



33%, 

which correspond approximately to the expected flavour tagging performance for this 
channel at the LHCb experiment 12 . In Table [2] we summarize the values of the physical 



parameters used to generate the toy data sets. 

We generate 500 data sets for different scenarios where we vary the values of the S- 
wave fraction Rs and its phase 5s and the weak phase —2j3 s . Each data set contains 
30000 signal events corresponding to approximately one quarter of a nominal LHCb year 
of2fb" 1 . 
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*Rs is fixed to when the S-wave component is neglected. 

Table 2: Values of the physical parameters used in the generation of signal decays and 
how these parameters are treated in the fit. 



We perform fits to each data set where 2f3 s , Rs, 5s, R\\, 5\\, R±, 5± are free parameters 
and all other parameters are kept fixed. We also perform fits where the S-wave component 
is present in the generated toy data, but ignored in the fit (Rs is set to 0) in order to 
investigate the bias in the determination of 2f3 s . 

The results of these fits for the statistical error and mean value of the weak phase — 2f3 s 
are summarized in Table [3j [4] and [5] for several different scenarios with — 2/3 s = —0.0368, 
— 2/3 s = —0.2 and — 2f3 s = —0.5, respectively. As an example, in Figure [T] we show the 
distributions of the fitted values of — 2j3 s for R s = 0.1, 5s = vr/2 and — 2f3 s = —0.5 for both 
the S-wave fraction Rs fixed to zero and Rs left free in the fits. In Figure [2] we show the 
distributions of the fitted values of Rs and the strong phase of the S-wave component 5s 
for the same case with Rs left free in the fit. It can be seen that when all parameters are 
fitted the results are unbiased, but when it is wrongly assumed that Rs = 0, the result 
for — 2(3 S acquires a bias with regard to the true input value. 

Figure [3] shows the bias in — 2f3 s from neglecting an S-wave component with Rs = 0.1 
and 5s = vr/2 versus the value of — 2/3 s used to generate the data sets. A linear dependence 
is observed, which demonstrates that the bias in — 2/3 s is proportional to the true value 
of —2j3 s . From Tables [3j [4] and [5] we observe biases for these scenarios which range 
from 7 — 17% in the measurement of 2/3 s if an S-wave component is present, but left 
unaccounted for in the fits. The bias moves the measured value of 2(5 S towards zero. 
This implies that the neglected CP-odd S-wave contribution has a bigger probability to 
be mis-identified as the CP-even longitudinal or parallel components than as the CP- 
odd perpendicular component. Therefore, although the bias from neglecting an S-wave 
contribution is unlikely to lead to false signal of new physics, it will cause a loss of 
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sensitivity to new physics. On the other hand, including the S-wave in the fit removes the 
bias in the central value of 2f3 s at a cost of an increase of less than 15% in the statistical 
error. 
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er(2/3 s ) = 0.045, Mean(2/3 S ) = 0.032 
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R s = 0.05, 6 S = tt/2 


<r(2/3 s ) = 0.048, Mean(2/3 S ) = 0.040 


cr(2/3 s ) = 0.045, Mean(2/3 S ) = 0.036 


R s = 0.05, 6 S = 


<r(2/3 s ) = 0.055, Mean(2/3 S ) = 0.038 


cr(2/3 s ) = 0.047, Mean(2/3 S ) = 0.032 



Table 3: Statistical errors and mean values of 2f3 s from 500 fits for different scenarios with 
2f3 s = 0.0368. The errors on a(2/3 8 ) and mean(2/3 s ) are approximately 0.003 and 0.002, 
respectively. The same data sets are used to obtain the results in the second and third 
columns. 
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Table 4: Statistical errors and mean values of 2(3 S from 500 fits for different scenarios 
with 2f3 s = 0.2. The errors on a{2f3 s ) and mean(2/3 s ) are approximately 0.003 and 0.002, 
respectively. The same data sets are used to obtain the results in the second and third 
columns. 





Float Rs in fit 


Fix R s to in fit 


Rs = 




<r(2/3 s ) = 0.051, Mean(2/3 S ) = 0.501 


Rs = 0.1, 6 S = tt/2 


cr(2/3 s ) = 0.059, Mean(2/3 S ) = 0.501 


cr(2/3 s ) = 0.053, Mean(2/3 S ) = 0.415 


R s = 0.1,6s = 
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R s = 0.05, 6 S = tt/2 
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Table 5: Statistical errors and mean values of 2(3 S from 500 fits for different scenarios 
with 2/3 s = 0.5. The errors on a(2/3 s ) and mean(2/3 s ) are approximately 0.003 and 0.002, 
respectively. The same data sets are used to obtain the results in the second and third 
columns. 
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Figure 1: Distributions of the fitted values of — 2f3 s for the scenario Rs = 0.1,5s = 
7r/2,2/9 s = 0.5. The left and right plots are obtained with or without fixing Rs to in 
fitting the data, respectively. 
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Figure 2: Distributions of the fitted values of Rs and 5s for the scenario Rs = 0.1, 5s = 
7r/2, 2f3 s = 0.5 without fixing R s to in fitting the data. 
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Figure 3: The bias in — 2f3 s from neglecting an S-wave component with R s = 0.1 and 
$s — 7r /2 versus the value of — 2f3 s used to generate the data sets. The bias is the difference 
of the mean of the fitted to the generated — 2f3 s values. A linear fit is superimposed on 
the graph. 
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4 Measuring cos2/? s 



In Eq. 2.1 and 2.2 one observes that the differential decay rates are invariant under the 



transformation 

(6\\ - 5 , 5± - S , S s - S , -2/3 B , Ar s ) «-> {S -5\\,ir + 5 - S ± , S -5 s ,ir- (-2/3 s ), -AT 

(4.1) 

As a consequence the measurement of 2f3 s is subject to a two-fold ambiguity, which is 
equivalent to cos2/3 s transforming into — cos2/3 s . A measurement of cos2/5 s including its 
sign would allow us to resolve this ambiguity. 

If the interference between the P-wave and S-wave amplitudes were to be significant 
in the 0(1020) mass region, we could use this effect to measure cos2/3 s , in the same way 



as BaBar measured cos 2(3 in B° — > J/0Kg7r° 13 . This requires measuring S,s — So, the 
strong phase difference between the S-wave and the longitudinal P-wave, as a function of 
the K + K~ mass in the 0(1020) mass region. When plotting this function, two branches 
are expected with each corresponding to a different solution for the weak phase (see 
Figure [1] left). It is straightforward to choose the physical solution since the phase of 
the P-wave Breit-Wigner amplitude is expected to rise rapidly through the 0(1020) mass 
region (dashed red curve in Figure [4] right), while the phase of the S-wave amplitude, 
which can be described either by a coupled channel Breit-Wigner function in case of an f 
contribution or by a constant term in case of a non-resonant contribution, is expected to 
vary relatively slowly (dotted green curve in Figure fright), resulting in 5s — So rapidly 
falling with increasing K + K~ mass (solid blue curves in Figure Eh. 
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Figure 4: An example to illustrate the dependence of the strong phase of the S-wave 6s , 
of the strong phase of the longitudinal P-wave So, and of their difference Ss — So, on the 
K + K~ mass. Left: the solid blue curve is the physical solution for Ss — So and the dashed 
black curve shows the mirror solution. Right: the dashed red, dotted green and solid blue 
curves are for S , Ss, and Ss — S , respectively. 



Below we use a Monte Carlo simulated toy data set to demonstrate the feasibility of 
this method in measuring the sign of cos2/3 s . We generate 30000 B° —>■ J/0K + K~ events 
in the K + K~ mass region between 1 and 1.05 GeV/c 2 , roughly corresponding to 0.5 fb -1 
of integrated luminosity. The P-wave and fo contributions are included coherently. The 
values of the parameters used to generate the toy data set are the same as in Table [2] 
except that we set — 2f3 s = —0.0368, and that the values of both Rs and Ss depend on 
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the K + K mass. The fo contribution accounts for about 10% of the total decay rate in 
the given mass region, as is shown in Figure [5] 




1005 1010 1015 1020 1025 1030 1035 1040 1045 1050 
K + K invariant mass (MeV) 



Figure 5: The data points correspond to the K + K~ mass distribution of a generated 
sample of B° — ► J/0K + K~ events including 10% fo contribution in the mass region. The 
dotted red curve indicates the fo contribution. 

The data sample is divided into bins in the K + K~ mass. For each bin i, two pa- 
rameters 5s,i and Rg,i are used to represent the average strong phase and the fraction of 
the fo contribution. Both sin 2/3 s and cos 2(3 S are treated as independent free parameters. 
Common free parameters sin2/3 s , cos2/3 s , R\\, R±, 5\\, 6±, T s and Ar s are used for all bins. 
Note that we still adopt the convention 6q = as only the relative phase differences in 
each bin can be measured. A combined fit to the time-dependent angular distributions of 
all the bins is performed to extract these free parameters. The fitted values of the strong 
phase difference 5s — $o versus the K + K~ mass are plotted in Figure [6} The two branches 
correspond to opposite values of cos2/3 s . Just as expected, the branch corresponding 
to the true solution decreases rapidly around the nominal 0(1020) mass. Choosing this 
branch leads to the unique solution 

sin 2(3 S = 0.043 ±0.05, cos 2/3 B = 1.05 ± 0.08 , (4.2) 
which gives the ambiguity-free result 

- 2(3 S = -0.043 ± 0.05 . (4.3) 

In this example, the measured — 2j3 s is separated from n — (— 2/? s ) by 13a, therefore 
the discrete ambiguity in 2(3 S is completely resolved. Although the actual measurement 
precision in cos 2/3 s will depend on the size of the fo contribution as well as background, 
the possibility to resolve the ambiguity in — 2/3 s using this method is very promising. 



5 Conclusions 

In the decay B° — > J/0K + K~ we expect that a K + K~ S-wave contribution in the narrow 
0(1020) mass region could be as large as 10%. The full differential decay rates for this 
decay including the S-wave contribution have been presented. We have considered a range 



S 




1000 1010 1020 1030 1040 1050 
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Figure 6: The fitted values of 5s — #o versus K + K~ mass are shown in red and black data 
points, corresponding to opposite values of cos2/3 s . The blue curve shows the dependence 
of 8s — So on K + K~ mass implemented in simulation. 

of scenarios which include S-wave components of 5% and 10%. We have shown that within 
these scenarios, if an S-wave component is ignored in the analysis, the measurement of 
the weak phase — 2/3 s would be biased by between 7% and 17% towards zero. We have 
demonstrated that by properly allowing for this S-wave component in the fit, an unbiased 
measurement of 2/3 s may be obtained with a slightly increased statistical error. Finally, 
we have shown that the interference between the K + K~ S-wave and P-wave amplitudes 
can be used to resolve the two-fold ambiguity in the measurement of the weak phase —2/3 s . 
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